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THERMAL NEUTRON ABSORBER BUILDUP IN IRRADIATED BERYLLIUM 

by Lawrence L Peters, Jr. 

Lewis Research Center 

SUMMARY 

The differential equations which describe the buildup of thermal neutron absorbers 
in beryllium have been solved. The solutions were used to calculate the change in the 
macroscopic absorption cross section of the Plum Brook Reactor as a result of the 
buildup of these poisons. This calculated value agrees within the experimental error of 
an actual observed change. The effects of variations in the factors on which the poison 
is dependent were studied and are presented in the form of graphs. 


INTRODUCTION 

When beryllium is exposed to fast neutrons, a chain of reactions takes place which 
results in the buildup of isotopes with high thermal neutron cross sections. For thermal 
reactors using beryllium as a reflector or moderator material these thermal neutron 
poisons can be significant, especially if the reactor must be shut down for several 
months or more. For example, after being shut down for approximately 100 days the 
PBR (Plum Brook Reactor) critical rod bank height at restart was about 0. 8 inch higher 
than normal. 

The buildup of poison in beryllium has been analyzed before by Morris (ref. 1) for 
continuous reactor operation. However, his solutions for the equations governing the 
poison buildup are not general because of the use of boundary conditions for continuous 
operation. 

In this work the differential equations governing the reactions were solved for more 
general boundary conditions. The change in the core absorption cross section was then 
calculated for cyclic reactor operation in the PBR. By using perturbation theory the 
corresponding change in reactivity was calculated. The analytical model was also used 
to study the effects of operation time, shutdown time, power, and fast to thermal flux 
ratio on the buildup of thermal neutron poisons in beryllium. 



SYMBOLS 


AK/K reactivity 

q q 

Nj number of Be nuclei per cm 

6 3 

Ng number of Li nuclei per cm 

3 3 

Ng number of H nuclei per cm 

3 3 

number of He nuclei per cm 

N? value of Ni at beginning of time interval 
52 change in macroscopic absorption cross section 

d 

3 -1 

A decay constant for H , sec 

Q 

<7 1 cross section for Be (n^, a) 

a 2 cross section for Li° (n g , a) 

O 

a ^ cross section for He (n g ,p) 

<p j fast neutron flvix > 0. 82 MeV, n/(cm^)(sec) 

2 

<p c thermal neutron flux, n/(cm )(sec) 

b 

cp* adjoint thermal neutron flux 


ANALYSIS 

Defining Equations 


The chain of reactions occurring in beryllium when in the presence of fast and 
thermal neutrons is as follows: 


Be 9 (n f , a) 


He 6 (j3")°- 84S 


Li 6 (n s ,.) 910b 


H fc 


(j3~)12. 5Y 

Z ' He 

(n g , p)5400b 


6 3 

Two of the isotopes formed, Li and He , have large thermal neutron absorption cross 
sections. To determine the concentration of these poisons, it was necessary to solve 
the set of simultaneous differential equations defining the reaction chain. These equa- 
tions are 
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= N lVl ' N 2^s a 2 


( 1 ) 


dN 

dt 


dNg 

dt 


N 2^s ff 2 


N 3 A+N 4^s a 4 


( 2 ) 


dN 4 

dt 


n 3 a 




•(3) 


The above equations do not consider loss of hydrogen or helium by diffusion; how- 
ever, for thick sections of beryllium (~1 in.) this is a negligible effect. The burnup 
is nearly zero, so for this analysis was assumed constant. The above equations 
were solved by using the Laplace transform and the complex inversion integral. The 
solutions give Ng, Ng, and as functions of time: 


2 Vz \ 1 


(4) 


(A+^s^Hj o ( N ? + N Vs CT 4 K Vl + N 2^s a 4^s a 2 

No = e <Ng + 

1 A + cp a 4 [X + <pAoa ~ (To)] (A- + cr 4 ) 


N °V2 


n 1<Ws<W4 


X * 'V°4 " n 2 } (X + <;> s a 4 ) 2 [x + f> s (a 4 - cr,)] 


s v 4 u 2' 


+ e 


-V2 1 N l^f CT l^s a 4 


N l^f CT l + N 2^s (a 4 - a 2 } | 


<P s o 2 [X + <P s (ct 4 - CTc,)] X + (P R (ct 4 - a 2 ) 


s v 4 u 2 


s v 4 2’ 


N 1 ^f g l < P s cr 4 
X + cp s o 4 


N% f a x + (No + N° + Njqa. 0 
— - — i - JLi - ^ViVt 

X + 


1 


1 


cp s a s ( ' X + A + V* 0 . 


s u 4 


s 4 


(5) 
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N 4 = e 


-(X + ^ s a 4 )t N 4 <p s a 4 


N 2 + N 3°)W 


A + < P s a 4 (A + <P s tf 4 )[X + <P s (ff 4 “ <t 2 )] 


n 3 °a 


N l^f a l^s a 2 A 


A + ^ g (a 4 - a 2 ) (x + + . ff? )] 


s' 4 " 2' 


, ^.^ Kvi N ol 

V2 


Nl < ^f a 1 At 


A + <P s (a 4 - a 2 ) A + <p c o 


s 4 


N 4 + N° + N° - nJ^otj 


^s a 2 X + <p « a 


s 4 


A 


A + <p„a 


s w 4 


( 6 ) 


With N 2 and N 4 determined, the increase in the macroscopic absorption cross section 
in the beryllium is given by 


6S a = N 2 a 2 +N 4 a 4 


( 7 ) 


A computer program was written to compute the poison buildup for cyclic reactor 
operation. (A Fortran listing of the code is given in the appendix of this report. ) For 
cyclic operation the values of N 2 , Ng, and N 4 at the end of an "on” or "off" portion 
of the cycle become the initial conditions for the next "off" or "on" portion of the cycle. 
During the reactor off time the values of <p s and <p f are negligible and taken as zero. 

To calculate the reactivity change corresponding to a given change in macroscopic 
cross sections, perturbation theory was used. (See for example, ref. 2.) The reac- 
tivity change is given by 

5S a / <^s dv ‘ 

AK/K = — (8) 

v I 2 f <p s dv 

The adjoint flux integral was obtained using the Twenty- Grand diffusion theory program 
(ref. 3). 
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Application to Plum Brook Reactor 


To determine the poison buildup in the PBR the beryllium reflector which surrounds 
the core on all four sides was broken up into seven regions. This was to account for the 
variation in flux around the reflector. Average fluxes for the seven regions were ob- 
tained from a two-dimensional neutron diffusion calculation. The value of cr. was deter- 
mined by flux- weighting the cross section over the interval 0. 82 MeV to 10 MeV. An 
average operation time and shutdown time per cycle was determined for the 5 years of 
operation prior to the 100 day shutdown. This information was input to the poison buildup 
computer program to solve for the poison atom densities. The reactivity corresponding 
to the change in absorption cross section was obtained from equation (8). The adjoint 
fluxes required to solve equation (8) were obtained from the neutron diffusion calculation. 
Table I gives some of the physical constants used for the PBR. 


TABLE I. - VALUES OF CONSTANTS USED 
FOR PBR CALCULATION 


N^, atoms /cm 3 

1. 12X10 23 

N®, atoms/cm 3 

0 

Ng, atoms/cm 3 

0 

N b , atoms/cm 3 

0 

a v b 

0.0214 

o 2 , b 

950 

°4’ b , 
sec 

5400 
1. 76x10® 

(p s , N/(cm 2 )(sec) 

<p t (>0. 82 MeV), N/(cm 2 )(sec) 

~5xl0 14 

~lxl0 14 

ton- da y s 

9. 45 

toff, days 

14. 95 


RESULTS AND DISCUSSION 

Comparison of Measured and Calculated Reactivity Change 

When the PBR was restarted following a 100-day shutdown, the critical rod bank 
height (RBH) was about 0. 8 inch higher than expected. In terms of reactivity this is 
about $1.35. The reactivity associated with the buildup of poisons in the beryllium re- 
flector was calculated to be worth $1.21. The estimated probable error in the predicted 


5 




or expected reactivity (RBH), and hence in the measured reactivity change, is ±30 percent. 
The estimated probable error in the calculated reactivity change is believed to be about 
20 percent. 


Effects of Variations in the Independent Variables 

Since there was a noticeable effect as a result of a 100-day shutdown, the question 
arose as to how serious the effect might be for longer shutdown periods. Figures 1 and 2 
show the poison buildup as a function of shutdown time following 5 years of cyclic opera- 
tion. For this period, which was the first 5 years of operation for the PBR, the average 
on time was about 10 days and the average down time was about 15 days. (Since that time 
this average has changed to about 20 days on and 5 days off.) 



Figure 1. - Poison buildup for first year of shutdown time 
for Plum Brook Reactor after 5 years of cyclic operation. 


6 




0 2 4 6 8 10 


Shutdown time, yr 

Figure 2. - Poison buildup as function of shutdown time 
for Plum Brook Reactor after 5 years of cyclic operation. 



Figure 3. - Poison buildup in Plum Brook Reactor beryllium 
for 20-day cycles. 


The present PBR core and fuel loading are such that the maximum excess reactivity 
in the control rods is about $15. 00. Thus from figure 2 it is seen that had the PBR been 
shut down for a 3 year period, it could not have been restarted without major core 
changes. 

Figure 3 shows the effect of operation time for 20-day cycles. The longer the on 
time per cycle, the more quickly the poison level rises. However, after about 1 year of 
operation the poison level for 5 days operation out of 20 and 15 days operation out of 20 is 
about the same. 

Figure 4 shows the poison buildup against operation time for several flux or power 
levels. (Figs. 4 to 6 are for continuous operation; however, the effects would be similar 
for cyclic operation. ) The higher the reactor power, the steeper the rate of poison 
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Figure 5. - Poison buildup as function of operation time 
with fast neutron flux ip f constant and thermal 
neutron flux <p s a parameter. 
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Figure 6. - Poison buildup as function of operation time with 
thermal neutron flux <p s constant and fast neutron flux 
(Of a parameter. 


buildup. However, after several years the poison level is independent of the power. 
After 10 years of operation the poison levels for the three different powers are within 
3 percent of each other. 

To examine the effect of the fast-to-thermal flux ratio, several combinations of the 
fast and thermal flux were calculated. The results are plotted in figures 5 and 6. Fig- 
ure 5 shows that, for a given fast flux, changing the thermal flux causes a change in the 
poison level and the rate of poison buildup. Figure 6 shows that, for a given thermal 
flux, increasing the fast flux by a factor of 10 causes an increase in poison by a factor 
of 10. 


CONCLUDING REMARKS 

Following a long shutdown of the PBR there was an unexpected change in the core 

6 3 

reactivity. This change can be accounted for by the buildup of Li° and He J in the beryl- 
lium reflector. 
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Although good agreement was obtained between the nominal measured and calculated 
reactivity change for the PBR, there was considerable uncertainty in both the measured 
value (because of uncertainty in the predicted critical rod position) and the calculated 
value. Chemical analyses are planned which will determine the atom densities of the 
reaction products in beryllium from the PBR core. This will give a second and more 
direct experimental check on the calculation of poison buildup in beryllium. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 20, 1970, 

120-27. 
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APPENDIX - IBM -1620 COMPUTER CODE FOR POISON BUILDUP IN 
BERYLLIUM FOR CYCLIC REACTOR OPERATION 


C POISON BUILDUP IN BERYLLIUM FOR CYCLIC REACTOR OPERATION 

5 PRINT 100 

100 FORMAT (57H1P01SON BUILDUP IN BERYLLIUM FOR CYCLIC REACTOR OPERATIO 
IN,/ ) 

READ 200, AN 1, AN2, AN3,AN4, PHI 1, PH I 4, NTS , TON, TOFF 
200 FORMAT (4E9.4, 2E9.2, I4,2F6.0) 

T=T0N*86400. 

DK= 1 . 76E-9 
SIGl=.0214E-24 
S IG2=950.E-24 
SIG4=5400.E-24 
PRINT 101, AN1 
PRINT 102, AN2 
PRINT 103, AN3 
PRINT 104, AN4 
PRINT 105 , PH I 1 
PRINT 106 , PH 1 4 


101 

FORMAT ( 23H 

IN I T I 

AL 

BE 

NO. DEN. 

IS, 

E10.2) 

102 

FORMAT ( 2 3H 

INITI 

AL 

LI 

NO. DEN. 

IS, 

E10.2) 

103 

FORMAT ( 23H 

IN I T I 

AL 

TR 

NO. DEN. 

IS, 

E10.2) 

104 

FORMAT ( 2 3H 

INITI 

AL 

HE 

NO • DEN • 

IS, 

El 0. 2/ ) 

105 

FORMAT (27H 

FAST 

(GRP. 

1) FLUX 

USED 

IS, E10.2 ) 

106 

FORMAT (27H 

SLOW 

(GRP. 

4) FLUX 

USED 

IS, E10. 2/ 


PRINT 117, TON 
PRINT 118, TOFF 

117 FORMAT ( 30H I RRADI AT I. ON TIME PER CYCLE IS,F6.2,5H DAYS) 

118 FORMAT (27H SHUTDOWN TIME PER CYCLE IS,F9.2,5H DAYS//) 

PHIF=PHI 1 

PHIS=PHI4 
1=0 
M = 0 

1 1=1+1 

A=AN1*PHIF*SIG1 
B =PH I S*S I G2 
C=PH I S*S I G4 
D=DK+C 
E = D-B 

EB=EXP ( -B*T ) 

ED=EXP(-D*T) 

BN2 = A/ B# ( 1 . -E B ) +AN2# EB 

BN3=ED*( AN3-( AN3+AN4)*C/D+( A+AN2*C)*B/E/D-AN2*B/E-A*B*C/D/D/E)+EB* 
1 ( AN2*B/E~ ( A+AN 2*0 /E + A*C /B/ E ) +A*C / D*T +( AN3+AN4 ) *C/D+ ( A+ AN2*C ) /D- A# 
2C*( l./B/D+l./D/D) 

BN4=ED*( AN4*C/D+( AN2+AN3 ) *B*DK/ D/ E-AN3*DK /E-A*B*DK/D/D/E )+EB* ( A/B- 
1 AN2)*DK/E+A*DK/D*T+( AN4+AN3+AN2- A* ( l./B+l./D) )*DK/D 
AN2=BN2 
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AN3=BN3 
AN4=BN4 
PN2=AN2*SIG2 
PN4 = AN4*S I G4 
DELS IG=PN2+PN4 
PRINT 1 08 t AN 1 
PRINT 109 » AN 2 
PRINT 110, AN3 
PRINT 1 1 1 * AN 4 
PRINT 1 1 2 t PN 2 
PRINT 113, PN4 
PRINT 114,DELSIG 
PRINT 115,1 

108 FORMAT ( 18H BE NO. DENSITY IS»E10.2) 

109 FORMAT ( 18H LI NU. DENSITY IS,E10.2) 

110 FORMAT ( 18H TR NO. DENSITY IS»E10.2) 

111 FORMAT ( 18H HE NU. DENSITY IS.E10.2/) 

112 FORMAT (25H DELTA SIGMA DUE TO LI I S , E 1 0 . 2 ) 

113 FORMAT ( 2 5 H DELTA SIGMA DUE TO HE IS,E10.2> 

114 FORMAT (25H THE TOTAL DELTA SIGMA IS,E10.2/> 

115 FORMAT ( 44H THE ABUVE VALUES ARE AFTER TIME STEP NUMBER , 14, // ) 
IF (M ) 2 , 2 , 3 

2 M = 1 

PH IF = 1 . E -20 
PHI S=l.E-20 
T=T0FF*86400 . 

I F ( I-NTS) 1,4,4 

3 M = 0 
PHIF=PHI1 
PH IS=PHI4 

T =T 0N*86400 . 

IF (I-NTS) 1,4,4 

4 PRINT 116 

116 FORMAT ( 19H END UF CALCULATION) 

GO TO 5 

END 
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